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Abstract 



in 

The spectral problem for 0{D) symmetric polynomial potentials al- 
P3 ■ lows for a partial algebraic solution after analytical continuation to nega- 

te ' tive even dimensions D. This fact is closely related to the disappearance 

^ , , of the factorial growth of large orders of the perturbation theory at neg- 

ative even D. As a consequence, certain quantities constructed from the 
perturbative coefficients exhibit fast inverse factorial convergence to the 
rvQ . asymptotic values in the limit of large orders. This quantum mechanical 

^ ' construction can be generalized to the case of quantum field theory. 

(N 

m 

O ■ 1 Introduction 

in 

f~^ ■ The 0(D) symmetric anharmonic oscillator 

o. 

o: ^ = ^E(p?+-?) + ^(E^O (1-1) 

i=l \i=l ) 

S^ . with a polynomial potential y(r^) is a good theoretical toy for testing various 

H \ ideas used in quantum field theory. In particular, one should mention: 

. 5^ ! 1) large orders of perturbation theory [DlllSllliniiailiaiS], 

2) fermion representations for boson theories in negative dimensions [TOlfTT] . 

3) hidden symmetries, quasi-exactly solvable models [TH|T3J|T3], correspon- 
dence between integrable quantum field theoretical models and ordinary differ- 
ential equations [151 HH |T71 [THl US] , correspondence between 0{D) symmetric 
and D = \ anharmonic oscillator systems HOI HT] |^ [S^ . 

4) Regge theory of complex angular momenta in nonrelativistic quantum 
mechanics [24 l [25 l [26] . 

There is a vast literature where the anharmonic oscillator is studied from 
various points of view. The aim of this paper is to bring some of these ideas 



together. The paper consists of two parts. The first part is devoted to the 
spectrum of the anharmonic oscillator at even negative D 

i? = 0,-2,-4,-6,... (1.2) 

understood in the sense of analytical continuation. A remarkable feature of these 
values of D is that a part of the spectrum can be found by solving algebraic 
equations. In principle, this fact is known since long ago [10, but as usual with 
exactly solvable models, one can arrive at the same results using many different 
ways, and the alternative derivations discussed here may deserve attention. 

The second part of the paper deals with the asymptotic behavior of large 
orders of the perturbation theory. As is well known, the perturbative series 
for the anharmonic oscillator is divergent because of the factorial growth of 
the coefficients. However, one can show that in even negative dimensions (jl.2p 
the factorial growth disappears and the perturbative series becomes convergent. 
Taken alone, this fact is quite natural and not too exciting. But this observation 
gives an interesting solution to another old problem. 

Let us return to the case of general D where we have the factorial growth 
of perturbative coefhcients. It is well known that the asymptotic behavior of 
large-fc orders of the perturbation theory is reached very slowly because of the 
0{k~^) corrections that are usually rather large in practically interesting cases. 
Therefore it would be rather interesting to construct quantities for which the 
asymptotic regime is reached faster than 0{k^^). This problem is especially 
acute in quantum field theoretical models where only several first orders of the 
perturbation theory can be usually computed. In this paper we construct and 
study special quantities whose approach to the asymptotic regime is factorially 
fast instead of the traditional slow 0(k~^) convergence to the asymptotic form. 
It is interesting that the same construction can be used in quantum field the- 
ory. In fact, this paper appeared as an attempt to understand some superfast 
convergence effects that are seen in quantum field theory [27 . 

Before starting with a detailed derivation of the exact solution for the spec- 
trum of the anharmonic oscillator at negative even dimensions D, it makes sense 
to describe the final result. For simplicity let us concentrate on levels with zero 
angular momentum / = 0. 

At negative even D a part of the discrete spectrum of energy E is given by 
roots of the algebraic equation 

RMiE)=0, (1.3) 

M--f = 0,1,2,3,... (1.4) 

where Rm{E) are certain polynomials of E depending on potential V . In this 
paper a rather elegant representation for these polynomials will be derived 

Rm{E) = det [J+ + V (2J_) - E] . (1.5) 

Here 

J±=Jl±lJ2 (1.6) 



are standard spin matrices for spin 

M D 1 3 , , 

j = - = -- = 0,-,l,-,... (1.7) 

Several comments should be made about the history of the exact solution 
for the anharmonic oscillator. The limit D —^ was considered by Dolgov and 
Popov [5] . From the point of view of representation (jl.Sp this case corresponds 
to spin j = with the polynomial Rq = — E leading according to eq. (jl.Sp to 
only level controlled by the exact solution E = Q. Although this case is trivial 
from the spectral point of view, the work of Dolgov and Popov [6| contains 
several interesting results. In particular, the wave function corresponding to 
the level E = Q was computed in Ref. [Ij. 

The general case of arbitrary negative even D was studied by Dunne and 
Halliday in Ref. [TU] where a method was suggested for a calculation of poly- 
nomials describing the spectrum. In principle, this solves the problem of the 
spectrum. However, some issues were not clarified. Testing some special cases, 
the authors of [lOj found an unexpected factorization of their polynomials into 
"elementary polynomials". In the framework of the current paper the "ele- 
mentary polynomials" of Dunne and Halliday are nothing else but polynomials 
i?Af(_E, f/) given by eq. ()1.5|1 . An important role in the disentanglement of the 
polynomial structure is played by a "hidden" sl{2) symmetry which stands be- 
hind the spin representation (jl.Sp and explains the miracles observed in Ref. 

m- 

In this paper we discuss several derivations of polynomials (II. 5p . One method 
described in section [2] uses an explicit expression for the Hamiltonian of the 
anharmonic oscillator in terms of differential operators obeying sl{2) algebra. 
Another way to polynomials (jl.Sp considered in Sec. 12.61 is based on the anal- 
ysis of a recursion relation for the coefficients of the power series for the wave 
function. In the third method (Sec. 14. 4p one performs analytical continuation 
of the I decomposition for the partition function. One more method (Sec. [5]) 
uses a fermion representation for the partition function continued analytically to 
negative dimensions. This approach is very close to the method of original work 
[TH] but we implement this method so that the sl{2) algebra remains explicit. 

An important aspect of the problem of the anharmonic oscillator in nega- 
tive dimensions is the precise definition of negative dimensions. It is interesting 
that the problem of analytical continuation to negative or complex dimensions 
is equivalent to the problem of analytical continuation to complex angular mo- 
menta which is a corner stone of Regge theory (Sec. 14. ip . 

As was already mentioned, at negative even D we meet interesting phenom- 
ena in large orders of the perturbation theory. The energy of the ground state 
of the anharmonic oscillator has a perturbative expansion 

oo 

£;(g,i?) = ^ii;W(i^)g'= (1.8) 

fc=0 

with coefficients E'^^^ (D) which are polynomials of D so that there are no prob- 
lems with analytical continuation of E'^^^D) to negative dimensions. In the 



case of the quartic anharmonic oscillator 

V{r^) = ^r^+gr^ (I.9) 

the large-fc asymptotic behavior of E^''^ (D) is described by formula [3] 

Ei^)iD) '=- i-ir^T [k + I) 3'=+^-^^ [1 + O(fc-i)] . (1.10) 

An interesting feature of this formula is that the RHS vanishes at even negative 
D because of the factor [T {D/2)] on the RHS (|1.10p . In section [6TT] we show 
that the factorial divergence of series ()1.8p disappears at negative even values 
of D so that the perturbative series has a nonzero convergence radius. 

Functions S^*^-* (D) are polynomials of D of degree k + 1 and they have fc + 1 
roots: 

E'-''\i^k,r) = (l<r<fc + l). (1.11) 

The last part of the paper is devoted to asymptotic properties of these roots in 
the limit of large k. One could wonder why we should care about the roots I'k.r- 
The answer is that these roots have a very interesting property: in the limit 
k ^ 00 some subsets of the roots are convergent to points D — —4, —6, —8, . . . 
and this convergence is factorially fast. In quantum mechanics this superfast 
factorial convergence may be an amazing but useless curiosity. However, in the 
context of quantum field theory analogous effects are important. Indeed, in 
quantum field theory one is usually limited to a rather small amount of pertur- 
bative terms. Four or five orders are usually considered as a great achievement. 
Therefore the existence of quantities with factorial convergence to the asymp- 
totic form instead of the slow 0{k^^) convergence of (ll.lOp is extremely inter- 
esting for field theoretical applications. The case of the anharmonic oscillator 
allows us to study this phenomenon in detail. 

2 0{D) symmetric Hamiltonians and s/(2) alge- 
bra 

2.1 Schrodinger equation 

The Z?-dimensional Hamiltonian with a spherically symmetric potential (|l.ip 
becomes in the polar coordinates 

In Z? > 2 dimensions, parameter I runs over nonnegative integer values 

/ = 0,1,2,3,... iD>2) (2.2) 



The degeneracy of the I levels is given by 

{2l + D-2)T{D + l-2) 



m{D,l) 



l\ 



V{D-1) 



(2.3) 



The cases D = 1, D = 2 are somewhat exceptional. But they still can be 
described by analytical continuation of (|2.3p in D (at fixed integer nonnegative 

For D = 2 and nonnegative integer Z, we obtain 



m{2,l) = 2~Sio, 

which can interpreted in terms of the ±1 degeneracy for I ^ 0. 
For D = 1, we have 

m{l,l) = ^(0 + Sn 



(2.4) 



(2.5) 



so that the only allowed values are / = 0, 1, and they correspond to the states 
with positive and negative parity respectively. 

After the separation of the factor r from the wave function one arrives at 



„-( 



1 



H-D = r-'Hy = -[- 



dr^ 



V- 1 d 
r dr 



+ V{r^) 



where 



V = D + 21. 
Thus the spectrum depends on D and I only via parameter T): 

Eni{D) = Eno{D + 21) = Eno{V) . 



(2.6) 
(2.7) 
(2.8) 



In the case Z? = 1 is the values ^ = 0, 1 correspond the states with positive and 
negative parity respectively. 
Introducing the variable 

^ 2 
we bring Hx> (|2.6p to the form 



(2.9) 



The corresponding Schrodinger equation is 



V d 



V{20 . 



-%-Ydc)+^(''^)-^"^^) 



^n(C)=0. 



(2.10) 



(2.11) 



Below we assume for simplicity that the potential V{2Q is polynomial in Q. 

Equation (j2.1ip was derived for positive integer dimensions T). However, 
we can use this equation for analytical continuation of eigenenergies EniT)) 
to arbitrary V. To this aim we solve equation (|2.11[) imposing the following 
boundary conditions on eigenfunctions iJjniC)'- 



1) ^n(C) niust be analytical at (" = 0: 

oo 

0(0 -Ep'^^"' (2.12) 

2) V'n(C) must decay at C — * +00. 

2.2 s/(2) representation for 0(D) symmetric Hamiltonians 

Let us define parameter 



J = -J (2.13) 



and operators 



To^-C^+J, (2.15) 

r_ = c. (2.16) 

These operators provide a representation of sl{2) algebra 

[T+,T_] = 2ro, [r+,ro] = -T+, [T_,ro] = r_ (2.17) 

with the Casimir operator 

C2 = ^ (r+T_ + T_T+) + ToTo , (2.18) 

[C2,r,]=0. (2.19) 

In representation (|2.14p - (|2.16p we have 

C2=j{j + l). (2.20) 

In terms of operators (|2.14p - (|2.16p Hamiltonian (|2.10p becomes 

Ht = T+ + V{2T^) . (2.21) 

In the case of the quartic anharmonic osciUator (II. 9|) we have 

ijQO = T++T^+ AgTl . (2.22) 

For integer or half-integer j one could expect a simplification of the problem 
of the anharmonic oscillator due to the presence of finite dimensional irreducible 
representations of sl{2). However, we meet a problem on this way: in the case 
of finite dimensional representations of sl(2) there exists a vector ip annihilated 
byT_: 

r_V' = 0, V't^O, (2.23) 



whereas for operator T_ (|2.16p equation 

T_V(C) = C^{z) = (2.24) 

leads to V' = [or to 'ip{z) — cS{z) if one allows for generalized functions]. In the 
next section we show how this problem can be circumvented in negative even 
dimensions. 

2.3 Effective spin Hamiltonian for even negative dimen- 
sions V 

Let us consider the case of positive integer or half-integer j 

J = 0,i,l,|... (2.25) 

According to eq. (|2.13p this corresponds to negative even values (|1.2p of param- 
eter T). These values are nonphysical. We assume analytical continuation in the 
sense of solutions of Schrodinger equation with boundary conditions described 
in Sec. O 

Let us consider the subspace Tij of functions /(C) obeying conditions 

/W(0) = (0<A:<2j). (2.26) 

This subspace is invariant under the action of operators Ta p.l4[) - (|2.16[) . 
Therefore we can reduce the action of operators Ta from the space H of differ- 
entiable functions to the factor space Ti-fTi-j ■ Expression (|2.2ip for Hamiltonian 
Ht shows that Tij is invariant also under the action of Ht- Therefore we can 
also reduce the action of Ht from H to H/T-Cj. The reduction to T-C/Hj is also 
possible for the spectral problem 

{Ht-E)^ = Q. (2.27) 

Certainly after this reduction only some finite part of the infinite spectrum 
can survive. Note that the factor space 'H/'Hj is isomorphic to the space of 
polynomials P2j[0 of degree 2j. The action of operators Ta (|2.14p - (|2.16p on 
these polynomials is described by the usual differentiation algebra extended by 
an additional formal rule 

^2j+fc ^ Q Jq^. ^^i (2.28) 

The action of operators Ta in the space of these polynomials is equivalent to the 
spin-j irreducible (2j + l)-dimensional representation of s/(2). This is obvious 
from the dimension 2j + 1 of this representation and from the eigenvalue j (j -I- 1) 
of the Casimir operator (|2.20p . Therefore after the reduction of the spectral 
problem (|2.27p to the factor space Ti/Tij we arrive at the matrix eigenvalue 
problem 

H.j = J+ + V{2J_), (2.29) 

[Hj -E)(j) = (2.30) 



where Ja are (2j + l)-dimensional matrices of the spin-j representation of sl{2). 
These sl{2) matrices Ja are connected by relations 



Jo = J3 , (2.31) 

J± = Ji ± iJ2 (2.32) 

with standard spin matrices Ji,2,3 for the spin-j representation of su{2). 
Thus the spectral problem reduces to solving the equation 

det {Hj -E) = 0. (2.33) 

Obviously 

R2j{E) = det (Hj ~E)^ det [J+ + y(2J_) - E] (2.34) 

is a polynomial of degree 2j + 1 so that we deal with an algebraic equation 

R2j{E)^Q (2.35) 

where j = — 1?/4 according to (|2.13p . 

One should keep in mind that matrix Hamiltonian (|2.29p describes only a 
part of the spectrum. Since sZ(2) matrices J± are real, polynomials R2j{E) are 
also real. Nevertheless the eigenvalues of Hj may be complex. 

In the case of quartic anharmonic oscillator (|1.9p matrix Hamiltonian (|2.29p 
reduces to 

Hj =J+ + J^+ AgJl . (2.36) 

2.4 Connection with quasi-exactly solvable problems 

The Hamiltonians of type (|2.36p appear in the context of quasi-exactly solvable 
(QES) potentials [121 [131 [14] . However, there is a certain difference between 
the sl{2) representation discussed here and in traditional QES problems. Our 
differential operators Ta (|2.14|) - p.l6p differ from the operators T^ used in QES 
problems: 

(2.37) 

(2.38) 

(2.39) 

Both sets of differential operators obey s/(2) algebra (|2.17p . In fact, T^ and T^ 
are connected (at least formally) by a Laplace transformation combined with the 
change of j — > —I — j. Nevertheless the two classes of problems are different. 
In particular, the methods discussed here allow us to compute a part of the 
spectrum but not the corresponding wave functions ip7i{C.) obeying equation 



r;: 


-^'^-^'t, 


n 




T'_. 


d 

~di- 



(|2.1ip . It should be emphasized that one has to distinguish two types of wave 
functions at negative V: 

1) solutions of Schrodinger equation (|2.1ip at negative even V, 

2) eigenvectors of the effective spin Hamiltonian (|2.29p . 

It is interesting that in the case of the sextic QES potential, Dunne and 
Halliday [TUj could trace the connection between the exact eigenf unctions in 
physical dimensions D = 1,2,3,... and the eigenvectors of the effective finite- 
dimensional matrix problem corresponding to _D == 0, —2, —4 . . . 

2.5 Example: harmonic oscillator 

In the case of the harmonic oscillator 

FHO(^2)^i^2^ (2.40) 

FH0(2C)=C, (2.41) 

the calculation of determinant ()2.30p is trivial: 



det ( J+ + J--E)= det (2 Ji - E) = det (2 J3 - E) ^ 2^^+^ JJ 



n=-j 



E 

n 

2 



so that solutions of eq. (|2.33p are 



(2.42) 



E = -2j,-2{j-l),...,2{j-l),2j 

V V v v , , , , 

= Y,Y + 2,...,-^-2,-- (I?<0). (2.43) 

It is instructive to compare this algebraic solution with the full spectrum of the 
harmonic oscillator. Labeling levels with n = 0, 1, 2, 3, . . . in each / sector, we 
can write 

EniiO,D) = ^ + i2n + l) . (2.44) 

Combining D and / into parameter V p.7p we arrive at the expression 

V 
E= — + 2n, 71 = 0,1,2,... (2.45) 

that has a trivial analytical continuation in P. For negative even V = — 4j 
the lowest 2j + 1 levels (|2.45l) obviously coincide with the solution of the spin 
problem ((2^ . 

In Fig. [T]we show the "Regge trajectories" (the relation between the ana- 
lytical continuation in D and Regge theory is discussed in Sec. 14. ip connecting 
physical values D = 1,2,3, .. . and with the special points D — 0, —2, —4, ... It 
is remarkable that any trajectory corresponding to a physical level after analyt- 
ical continuation to negative D sooner or later enters into the domain described 
by the roots of polynomial (|2.42p . 




Figure 1: Trajectories continuing the levels of the harmonic oscillator analyti- 
cally in T) from physical values V — 1,2,3,... (marked with crosses) to arbitrary 
values of V. The levels at I? = 0, —2, —4, . . . that are described by roots of poly- 
nomials i?2j (j = —D/2) are marked with circles. 



2.6 Power series and recursion relations 

The analytical continuation of levels En{'D) to negative V uses solutions 0(C) 
(|2.12p of Schrodinger equation (|2.1ip that are analytical at ^ = 0- For polyno- 
mial potentials 



obeying condition 
let us define 

Then 



V{20 = Y.WkC\ 



Wk = if A: < or k > L 



V{0) = 

Wk if i < k < L, 
Uk ^ { -E if fc = 0, 
otherwise . 



V{20-E= J2 "feC 



(2.46) 

(2.47) 

(2.48) 
(2.49) 

(2.50) 
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Inserting this series and expansion (|2.12p into cq. (|2.1ip one obtains the recur- 
sion relation for the coefhcients pk of series (I2.12p : 



L 

ak+l{T>)pk+l = X! '^mPk-m , (2-51) 

rn— 

ak{V) = fc ( fc - 1 + -v\ . (2.52) 

If 13 7^ 0, —2, —4, . . . then we can solve this equation iteratively starting from 
some po 7^ and obtain 

Pk+i = -pf-. (2.53) 

where Sk+i ({um},I?) are polynomials of Um- Once the series (I2.12p is con- 
structed, one has impose the boundary condition on 0(C) at C ^ +oo and this 
will fix the spectrum. 

The case V = — 4j = 0, —2, —4, ... is exceptional because in this case at 
k = —'D/2 = 2j we have a zero in the denominator on the RHS of (|2.53p . 
Therefore instead of the determination of p2j we will arrive at the condition 

S2j+i({wm},-4j) = 0. (2.54) 

Since S2j+i ({um},2?) is a nonzero polynomial in E (via uq = —E), this condi- 
tion does not hold for the most of energies. In fact, in the case V = 0, —2, —4, . . . 
the majority of eigenfunctions have power expansion (j2.12p with 

Po = Pi • • • = P2J+1 = 0, P2J+2 ¥" . (2.55) 

However, if condition (|2.54p holds then we can proceed with iterations and 
construct the series (|2.12p . Thus in the special case ()2.54p we have two solutions 
p.l2p of Schrodinger equation (|2.1ip analytical at C = 0: one solution with 
Pq y^ 1 and the second solution (|2.55p . Thus all solutions of Schrodinger equation 
(|2.12p are regular in the case (|2.54p . This means that taking an appropriate 
linear combination we can satisfy boundary conditions at infinity so that values 
of E obeying equation (|2.54p automatically belong to the spectrum. 

It is easy to see that this construction is equivalent to the spin problem 
(|2.33p . Indeed, recursion relation (|2.5ip is nothing else but vector equation 
(|2.33p written in the basis diagonalizing matrix Jq. Now we understand that 
polynomials S2j"-i-i {{wm}, E,—Aj) appearing in eq. (|2.54p must coincide up to 
a constant with polynomials R2j {E) defined by eq. (|2.34p 

S2J+1 ({um}, -4j) = c, det [J+ + V (2J_) - E] (2.56) 

with some coefficients Cj. 

Let us derive equation (|2.56p more carefully and compute coefficients Cj. 
First we define a (TV -I- 1) x (A^ + 1) matrix 

Ci'^,'' ^ Uk-n ~ ak5k^n-i l<k,n<N+l. (2.57) 
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According to eq. (|2.49p we assume that u^: = if fc < 0. 
Then recursion relation (|2.5ip takes the form 



Af+l 



Y,ci^^Pn-i^6k,N+iaN+iPN+i {l<k<N + l). (2. 



Matrix (|2.57p is quasitriangular. Its determinant can be computed by making 
with its rows the same hnear manipulations as in the recursive calculation of 
Sk+i ({um},I?) in eq. (|2.53p . As a result, one obtains 



JV 



det 



C(^) ({«„,}, V) = SN+i ({«„,}, V) Yl ak{V) . 



(2.59) 



k=l 



Now we set V = -4j, N = 2j. Using eq. (|232|) . we find 



2j 



n«fc(-2j-) = (-l)'M(2j)!f 



fe=i 



Thus 



det 



C(2^) {{u^}, -Aj)] ^ (-1)2^- [(2j)!] S2,+i ({«„}, -4j) 



Therefore algebraic equation (|2.54p takes the form 



det 



C(2j")(K„},-4j) 



= 0. 



(2.60) 



(2.61) 



(2.62) 



In practical calculations one can compute detC^^-^-' either directly using defini- 
tion ((^37)1 

detC(2j)- det [wk^n-k{k-l-2j)Sk.n-i-ESk„] (2.63) 

l<k,n<2j + l 

or solving recursion relations iteratively and finding in this way S2j+i {{um}, — 4j). 
Both methods lead to the same results: 



det C^" 
detC^i 
det C(2 
det C(3 
detC^-^) = -£""• + 20E'' 



= -E, 

= -E^ + 4?«2 + 4Ewi , 

= £;■*- lOE^wi + 9wl ~ 24w2E - 86103 , 



(2.64) 
(2.65) 
(2.66) 
(2.67) 



wi - 64£;wf + 8Aw2E^ - 192u;iW2 + 288Ew3 + 576^4 . 

(2.68) 



Note that recursion relation (|2.58p for N = 2j is nothing else but the equa- 



tion 



[T+ + V{2T^)-E]^iC)=0 



(2.69) 
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written in the basis ('', k = 0,1,2, .. . 2j. This basis coincides with the standard 
basis for spin-j representation up to the subtleties of the normahzation and 
enumeration. Obviously the determinant of matrix C^'^^' is invariant with the 
respect to this trivial change of the basis. Therefore 



detC^^^') =det[J+ + F(2J_) 



^1 



(2.70) 



where the spin j representation is assumed on the RHS. Combining this result 
with eq. (|2.62[) we see that we have obtained old equation (|2.29p . 

Certainly this "new derivation" of eq. (|2.29p is nothing else but an explicit 
detailed version of the compact arguments of Sec. 12.31 presented there in terms 
of the sl{2) algebra acting in the factor space Ti/Tij. 



2.7 Properties of polynomials R2j{E) 

Comparing eqs. (|2.34p . (|2.63p . (|2.70p we see that we have three equivalent 
representations for polynomials i?2j {E) 



R2j{E)=dct[J+ 



det 

l<fe,ri<2j + l 



y(2j_) -E] 

[Wk^n - fc (fc - 1 - 2j) dk, 



-1 - E6kr, 



= i-ir^i2j)\]s2,+ii{u^},-4j). 



(2.71) 
(2.72) 

(2.73) 



Representation (j2.7ip is written in terms of spin-j matrices J± = Ji ± i J2 and 
explicitly expresses the sl{2) symmetry standing behind this construction. Rep- 
resentation (|2.72p written in terms of coefficients Wk of the polynomial potential 
V (|2.46p is nothing else but eq. (|2.7ip transformed to a basis with a nonstandard 
normalization which allows us to get rid of cumbersome square roots appearing 
in the traditional expressions for the matrix elements of J±. Representation 
(I2.73P is based on the iterative solution ()2.53p of recursion relations (12.511) and 
is useful for computer calculations of polynomials R2j{E). 

The higher coefficients of polynomials R2j{E) can be computed using the 
large- i? expansion: 



det [E- J+- y(2J_)] = E^^+^ exp <^ Trln 



J+ + T/(2J_) 



E 



^2J+i cxp I - ^ -E-"Tv [J+ + y(2J_)]" I . 



(2.74) 



We have 



Tr[J+ + F(2J_)] =0, 
Tr|[J+ + T/(2J_)]H ^ 2Tr[J+F(2J_)] = 2wiTr (J+J_) , 

Tr [J+ + V{2J^)f =: 3Tr \{J+f F(2J_)j = 3w2Tr \{J+f (J_)^j (2.75) 
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where Wk are coefficients of the polynomial potential (|2.46p . A straightforward 
spin algebra gives 



Tr 
Finally we obtain 



Tr (J+ J_) = -jij + l)(2j + 1) , (2.76) 

{J+f {J-f] - ^3{j + l)(2j + l)(2j - l)(2j + 3) . (2.77) 

10 



- ^j(j + l)(2j + l)(2j - l)(2j + 3)W2E'^-^ + 0{E'^-'), 

(2.78) 
In the case of harmonic oscillator (|2.40p we have according to (|2.42p 



Rfp{E) = det (J+ + J--E) = 22J+1 Y[ 



E 

n 

2 



(2.79) 



2.8 Symmetry D ^A- D 

Relation 






shows that equation (|2.11|) has a symmetry 

X>^ 4-1? . 



)(C) 



(2.80) 



(2.81) 



(2.82) 

However, one should be careful about the boundary conditions. Indeed, for 
V < 2 the factor C~i+(^/2) is singular so that regular solutions 4>t>{C,) may 
correspond to singular solutions of (/)4_x)(C)- In other words, some part of the 
T) <2 spectrum may be lost when one turns from I?<2to4 — I?>2. 

According to results of Sec. 12.61 ai V — 2M = 0, —2, —4, . . . there are two 
types of regular solutions (|2.12p 

1) solutions with pa ^ 0, 

2) solutions with po ~ . . . ^ pM ~ 0, pm+\ t^ 0. 

Obviously solutions 0-d(C) of the first type with po 7^ generate singular 
functions (/)4_-p(C) (|2.8ip . Therefore the spectrum of the 4 — 1? problem does 
not contain that part of the T) = —0, —2, —4 spectrum which corresponds to 
wave functions (|2.12p with po = 0. We see that under the change I? — > 4 — I? 
we lose exactly those states which are described by the roots of polynomials 
Rm{E'). 

This symmetry I? ^ 4 — I? is illustrated in Fig. [2] It should be stressed that 
the symmetry I? —> 4 — I? is relevant only for integer even values of 2?. 
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3 Quartic anharmonic oscillator 

3.1 Polynomials R2j{E,g) 

In the case of the quartic anharmonic oscillator (|1.9p we have 



V{20 = (: + Age^w,+W2C\ (3.1) 

u;i = l, W2^4g. (3.2) 

Now one can use one of representations (|2.7ip - (|2.73p and compute i?2j (E) . In 
order to avoid large coefficients it is convenient to define 

R2j{E,g) - (-2)-2J R^^{E) = (-2)-2J det (j+ + J_ + igJl - E) . (3.3) 



At = we have according to (]2.42p 






fc=0 ^ ^ 




(3.4) 


Let us define 






R'^^{E,g) = R2,{E,g) - R2j{E,{)) . 




(3.5) 


The first polynomials i?2j {E, g) are 






R[{E,g)={), 




(3.6) 


R'2{E,g)^Ag, 




(3.7) 


R'^{E,g)^l2Eg, 




(3.8) 


R'^{E,g) ^3{^16 + 7 E^)g, 




(3.9) 


R'^{E,g) = 4(7^3 - 55£: - 2005)g , 




(3.10) 


R'&{E,g) = ^(7i?4 - 124i?2 + 192 - imQEg)g , 




(3.11) 


R'j{E,g) = ^(7£;^ - 230£;3 + 1183S ~ ZQb&E'^g + 


109765)3 . 


(3.12) 


Now we find using eas. (|3.4p and (|3.5p 






i?o(^,5) = -^, 




(3.13) 


Ri{E,g)^-\{E''-\), 




(3.14) 


i?2(i?,.g) = -Ji^(i?'-4) + 45, 




(3.15) 


i?3(£;,.9) = -i(i?' - 9)(ii;2 - 1) + 125i?. 


(3.16) 
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3.2 Ground state at D = -4 

The ground state has Z = so that the case of the ground state in D = —4 
dimensions corresponds to V = D = —4. The spectrum is described by the 
roots of polynomial R2{E,g) (|3.15p 

R2iE,g) = (3.17) 



I.e. 



Let us define 



E^ - 4:E - 16g = . 



E=^s, (3.18) 

h = 3^^^g. (3.19) 



Then 



e^~3e-2h = 0. (3.20) 

The solutions can be found using Cardano formula: 

e= \/h+^i^^+^h-^/h^~l. (3.21) 

Keeping in mind applications to the case of the perturbation theory in small g 
it is convenient to transform this expression to the form 



-V3 + i' 



2Re 



y/l-h^- ih 



(3.22) 



which allows for the perturbative expansion 

£ = -V3+J + ... (3.23) 

o 

This leads to the perturbative expansion for the energy of the ground state in 
D = —4 dimensions 

Eig, -4) = -2 + 2g + 3c/ + Sg^ + 1^/ ^ gg^s + ^g^ + 1536/ + . . . (3.24) 

4 ■ lb ■ 

In Sec. [6] we will compute the asymptotic behavior (|6.23p of this series. 

4 Analytical continuation to negative D 

4.1 Negative V and Regge theory 

We have already made some comments about the role of Schrodinger equation 
(j2.11|) for the analytical continuation in V. Remember that parameter T> ()2.7p 
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is built of D and I so that the problem of analytical continuation in V has two 
equivalent formulations: 

1) analytical continuation in D at fixed /, 

2) analytical continuation in I at fixed D. 

The second approach allows us to use some results from Regge theory of 
complex angular momenta [231 US [5nj . Although the traditional Regge theory 
deals with the scattering problem whereas we are interested in polynomial po- 
tentials, the analysis of the behavior of the solutions of Schrodinger equation 
(|2.1ip at C ^ is the same in Regge theory of the potential scattering and in 
our case. In particular, in Regge theory for potentials regular at r = 0, the 
Regge poles are at points 

3 5 7 ,^^^ 

' = -2'-2'-2--- (4-^) 

Combining this with the general formula (|2.7p applied to the case Z? = 3 

V = 3 + 2l {D^3) (4.2) 

we see that ^-points ()4.1|) correspond to negative even values V = 0, —2, —4, —6, . . 

From the point of view of the historical perspective it is also interesting 
that analytical continuation in I for potentials growing at infinity, e.g. for the 
oscillator energies (j2.44p , was discussed in the context of naive models for quark 
confinement p51l26j . 

Although the discussion of analytical continuation in terms of / is preferable 
from the point of view of explicit connections with Regge theory, we choose 
the D representation. The D language is natural for the analysis of analytical 
properties of the partition function (Sec. 14. 4[) and for the fermion representation 
at negative D ^. 

4.2 D dependence of levels 

We have already discussed the dependence of energy En {T>) on T) for the har- 
monic oscillator (see Fig. [1]). If we switch on a small anharmonicity then the 
trajectories will be deformed but qualitative features will be preserved. How- 
ever, in the case of a strong deviation from the harmonic regime (or at larger 
values of |I?|) some new phenomena may occur. In Fig. [2] we show the T) de- 
pendence of levels for the quartic anharmonic oscillator ()1.9p with ,g = 1. We 
see that the trajectories of the ground state and of the first excited state meet 
together at I? « —2.6 . Polynomials i?2j have 2j + 1 roots but some of these 
roots may be complex. In Fig. [2] we see that this happens for polynomials i?2, 

i?3, Ri. 

4.3 Various aspects of the problem of analytical continu- 
ation in V 

The problem of analytical continuation in T) can be studied for various quanti- 
ties: 
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Figure 2: Trajectories continuing the levels of the quartic anharmonic oscillator 
with 5 = 1 analytically in T) from physical values P = 1, 2, 3, . . . (marked with 
crosses) to arbitrary values of T). The levels at I? = 0, —2, —4, . . . that are 
described by roots of polynomials i?2j {] — —D/2) are marked with circles. 
Note that the trajectories of the ground state and of the first excited state meet 
together at D « —2.6 . Some polynomials R2J have complex roots so that the 
number of real levels associated with these polynomials is smaller than 2j/ + 1 . 
The dashed horizontal lines connect levels related by the symmetry D ^ A — D 
(for even D). 
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1) analytical continuation of separate energy levels 

E,,iiD) = E,,o{D + 21) = E,,o{V) , (4.3) 

2) analytical continuation of coefficients En (V) of the perturbative expan- 
sion 

oo 

£;„o(2?,5) = E^io^(^)3' (4-4) 

fc=i 

for anharmonic oscillator, 

3) analytical continuation of the partition function 

Z{P, D) = Tr exp {-(3H) . (4.5) 

Let us comment briefly on the analytical properties of these quantities. 

The situation with analytical continuation of the perturbative coefficients 
i?„g (I?) in expansion (|4.4[) is simple because -E„q (I?) are polynomials of V. 
However, the asymptotic nature of the perturbative series (|4.4p does not allow 
us to draw simple conclusions about the analyticity in V of exact levels Eno (T>, g) 
from the trivial analyticity of _E„q (V) . 

The simplest example of analytical continuation in D provides the harmonic 
oscillator. Its partition function trivially factorizes into D components 

zW(/3,i?,0)=[zW(/3,l,0)]'' (4.6) 

with an obvious analyticity in D. The case of anharmonic potentials will be 
considered in the next section. 

4.4 Analytical continuation of the partition function 

Now let us study analytical continuation of the partition function 

Z{l3,D) = Tvcxp{-l3H) (£> = 1,2,3, .. .) (4.7) 

in D. At physical values oi D — 1,2,3,4,... one can use the I decomposition 

oo 

Z(/3, D)^Y. ™(^' exp [~l3Eni{D)] . (4.8) 

Here m{D, I) is the degeneracy of the Z-levels in the £)-dimensional space given 
by eq. ([2?3ll . Using eq. ([2?8ll . we find from eq. (|48ll 

Zip, D)^J2 "^(-O' ^) I J2 exp [-(3EMD + 21)] I . (4.9) 

1=0 [ n J 
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Let us define 

oo 

z(/?, D) = Y, cxp [-(3EMD)] . (4.10) 



n=0 

Then 

OO 

Z(/3,i:») =^m(i:>,02(/3,£) + 20. (4.11) 

We want to use this series for analytical continuation of Z{/3, D) to negative D. 
It is easy to check that for positive integer / = 0, 1, 2, 3 . . ., function ■m{D, I) is 
analytical in D in the whole complex plane. Analytical continuation of m{D, I) 
to D = 1 and D = 2 was already discussed in Sec. 12.11 Now let us consider the 
general case. Using the relation 

r(z)r(i - z) = ^^^ , (4.12) 

siuTT^; 
we can transform eq. (|2.3p to the form 

, „ , ,,(21 + D -2) r(2-D) 

Combining representations (|2.3p and (|4.13p . we see that m{D, I) is regular for all 
D. The regularity of m{D, I) allows us to use representation (|4.1ip for analytical 
continuation of Z{D) to negative D. 

At negative even values D — ~2M — 0, —2, —4, ... we have an important 
simplification. Indeed, according to eq. (|4.13p 

^(-2M,0 = (-l)'+^(^^-^f-^^ r(2 + 2M) 

^ ' ^ V y ^1 V(Z + 2M~l) ^ ' 

We see that 

m{-2M,l) = if />2M + 3. (4.15) 

Therefore series (I4.1ip reduces to a finite sum: 



(2 + 2M-n! 
1=0 \ < / 



(4.16) 



Now we change the summation variable from I to 

k = l-M -1. (4.17) 

Then 

7(f^ OAn \^' (-l)^-^'^2fc[(l + 2M)!] .._„,. .,.„. 
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Note that k — does not contribute to this sum. Combining the contributions 
of k and — fc we arrive at 

Z(3 2M)-'^' (-l)^-^^2fc[(l + 2Af)!] 

X [z(/3,2 + 2fc)-z(/3,2-2fc)] . (4.19) 

Now we use the D ^ A — D symmetry (|2.8ip connecting the levels of z{P,2 + 2k) 
and z{P, 2 — 2k). As was explained in Sec. 12.81 this symmetry is partiah when 
one passes from Z) = 2 — 2fc to I? = 2 + 2A: one loses the states described by the 
solutions of the algebraic equation Rk-i{E) = 0. Therefore 

z(/3,2-2fc) = z(A2 + 2fc)+ [Trexp(-/3iJ,^_(fe_i)/2)] (4.20) 

where Hj is the matrix Hamiltonian (|2.29[) corresponding to the spin 

fc-1 



(4.21) 



Now eq. (|4.19p takes the form 



Z(B 2M) V (-l)'--"^+^2fc[(l + 2M)!] 

Z{P, -2M) - 2^^ ^^^ ^ ^ ^^,^^ ^ ^^ ^ ^^, (4.22) 

xTrexp(-/3iJ,^(fc_i)/2) . (4.23) 

Let us replace the summation variable k to j (|4.2ip 

Z{(3,-2M) = y (-l)2^-^^-^^^l±Mi±Mi-Trexp(-/3iJ,) . (4.24) 
^^' ' ^^ ' {M - 2j)l{2 + M + 2j)l i'y >^ J' V ; 

Here the summation runs over integer and half-integer j. 

5 Fermion representation 

5.1 Derivation 

Let us start from the case of the quartic anharmonic oscillator (|1.9p . We can 
write the path integral for the partition fmiction 



Z^^{P,D,g)^ I Dxexpl- I dt 



<2(0)=<2(/3) 



1=1 \j=i / 



(5.1) 
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One can use the standard trick with the auxihary field a 



Z'^''{l3,D,g)^ I Dx I Da 

q{0)=q{T) 

exp<-/ dt -^Xi(-dt +l + i^ 



Da 



Det 



(-* 



+ l + i 



-D/2 



exp 



, a 
dt — 
2 



For negative even D 

D = -2M 

we can use the fermion representation 

-D/2 



(5.2) 
(5.3) 



Det 



(-«■ 

Dtp^ Dtp exp 



l + iW8qa 



l3 



dt 



M 



Y,i^t{~d^ + ^ + ^V^9'^)^' 



(5.4) 



Integrating out the field a, we obtain 
Z<i^{P,D,g)''-^'^' [DtP+DtP 

r0 



X exp 



dt 



M 



M 



J2 [{dti't) iM^) + i't^^] + 45 J2^t^^ 



(5.5) 



The same trick can be repeated (at least formally) for an arbitrary polyno- 
mial central potential V{r'^) if one uses more intricate integral representations 
involving several auxiliary boson fields a. In this way one arrives at the fermion 
representation for partition function (|4.5p with an arbitrary polynomial V: 



ZiP^D) ""-=''' J 



DTp+Dip 



X exp ■ 



P 



dt 



M 



M 



Y,[{dttPt) {^t^^)] +V [2Y,^t^^ 



(5.6) 



Now we introduce the auxiliary fermion field Xi'- 

J Dx+Dx exp ll dtJ2 [{dti^t) X^ " xt i^tiJ^) - xtx^] \ 



exp 



dtY^{dti^t){dti^i) 



(5.7) 
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Then 



Z{P,D)''^=^^ / D4>^DijDx^Dx 



X cxp <^ I dt 



M 



i=l 



Next we change the notation for the integration variables: 

ipi -> an , 



Vi - 


-" ^a ' 


Xi ~ 


-> ai2 , 


xt- 


-^4 



(5.8) 



(5.9) 
(5.10) 
(5.11) 
(5.12) 



so that 



Z{p, D) ^^=2^ / Da+Da 



X exp ■ 



/3 
/O 
M 2 



(ft 



M 2 



H Y^ ""i^ ( '^"'3^* + 9 (^1 + *'^2)a/3 ) «^/3 



«=1 Q,/9=l 



-^ I XI XI "i^ ('^1 ^ *'^2)„/3 a*/3 
i=l a,/3=l 



(5.13) 



where Ca are standard Pauli matrices. 

The integral on the RHS corresponds to the effective fermion Hamiltonian 

M 2 -. I M 2 \ 

i=\ a, 13=1 V=^ aJ3=l J 



with anticonimutation relations 

{flma, an/3} = |a+„, a+^l == . 



Thus we have derived the formula 

D^-2M 



Z{P,D)''^=^Tt [e-'^^-Pf] 



(5.14) 

(5.15) 
(5.16) 

(5.17) 



which should be understood in the sense of analytical continuation of the bosonic 
partition function Tr e^^^ from positive integer dimensions to negative even 
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dimensions. On the RHS we have a trace in a 2^*^ dimensional vector space 
Pp is the operator counting the fermion parity of the state: 



Nf=Y, "mO*' 



(5.18) 
(5.19) 



Pp appears in eq. (|5.17p because the integral over Grassmann variables (|5.13p 
inherits periodic boundary conditions of the original boson integral (|5.ip . 
Now we define 



M 2 



1=1 a,/3=l 



^'^ " 9 XI Yl ^"^ ('^a)a/3 a^0 

Obviously 

Introducing sl{2) generators 



Ta,Th 



— ^^ahc-^ c 



we can write Hamiltonian (|5.14p as 

This effective Hamiltonian has a form similar to form (|2.2ip . 
5.2 Degeneracy factors 



(5.20) 

(5.21) 

(5.22) 
(5.23) 

(5.24) 



Note that in the first derivation of the spin Hamiltonian Hj (|2.29p via the sl{2) 
expression Ht (|2.21[) we worked with the "effective dimension" V given by (|2.7p 
keeping in mind that the spectrum depends on D and I only via V. On the 
contrary, in the derivation based on the fermion Hamiltonian Hp ()5.24|) we can 
trace the contributions of levels with different I. Strictly speaking, Hp (I5.24|) is 
not completely equivalent to Hj (|2.29p because Hp acts in the 2^*^ dimensional 
space corresponding to the reducible representation of sl{2) (or su{2) if one 
prefers a more familiar physical interpretation) 



M 



(g) [o]e[o]e 



(5.25) 



Here we use notation [j] for irreducible representations corresponding to spin 

j. Each factor [0]e[0]e [i] 

sector: 

'1" 

|0> 



in (I5.25P is associated with the states in the i-th 

«ilo) 



[0] 
«a«.+2|0) ^ [0] 



(5.26) 

(5.27) 



(5.28) 
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Tensor product (|5.25p can be decomposed in irreducible representations [j] 

M , r , -1 \ J^-f/2 

2 



/ rii\ ' 

(2)[0]e[0]e - ]=@n{j,M)\j] (5.29) 



i=l 

with degeneracies 



3=0 



, _, 2 (2? + 1) T(2M + 2) , , 

^(^' ^^) - r(M-2, + i)r(3 + M + 2,) ^'-''^ 

which can be computed using the recursion relation 

n(i, Af + 1) = 2n(j, M) +n(j- ]^,m\ +n(j + i, ilA . (5.31) 

Now we apply decomposition (|5.29p to the calculation of the partition func- 
tion (|5T7)l : 

M/2 

Z(/3, i?) ^^=''' Y. {-^?'-''<J^ M) [Trcxp (-/3iJ,)] . (5.32) 

3=0 

Inserting n{j, M) from eq. (|5.30p . we arrive at the result coinciding with (|4.24p . 
Thus fermion representation (|5.6p for the partition function leads to the same 
result as the direct calculation of Sec. 



5.3 Connection with the results of Dunne and Halhday 

If one uses representation 



(5.33) 
(5.34) 
(5.35) 
(5.36) 



2—1 \% — l / 

M 2 / M 2 

^F=^ Yl 2°'tc (f^i + «'72)„^ a,^ + ^ XI Y^ 4 ('^i - ^^2)„^ a.p 

i=l Q,/3=l \4=1 a, 13=1 









4 = 


^2J-l , 








0-12 = 


d 




afe-i 








Oil = 


^2J , 








4- 


9 

^92^ 


Hamiltonian 


Hf 


takes the form 






Hr^ 


\ 


d d 


— + V 



(5.38) 
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^^-- 


- ^2'i-1^2i , 


Ci^k-- 


= -&6 , 


m' = 


= 



used by Dunne and Halliday [TU] ■ The advantage of this form is the possibility 
to work with the pseudo-boson variables 

(5.39) 
(5.40) 
(5.41) 

and to look for the eigenstates of Hp in the form of polynomials ^({^i}) ^^ ^i- 
In our representation these polynomial wave functions correspond to the 
states of the form 

P({a+aa})|^) (5.42) 

where P are polynomials of a^Oii, and \n) is fixed by the conditions 

a+\n)=ai2\^)=0 {l<i<M). (5.43) 

This corresponds to working in the subspace of states obeying the constraint 

2 

Y, a+^a,f3i> = i> (1 < i < M) . (5.44) 

/3=1 

In terms of decomposition (|5.29p this constraint selects states (j5.26p in each 
«-sector. Therefore states (I5.42p belong to the subspace: 

M r -, -1 *V2 

fi,{j,M)[j\, (5.45) 

] = {M/2} 

where {a;} stands for the fractional part of x, and the degeneracy factors 

follow from the recursion relation 

n(j, Af + 1) = n ("j - 1, mJ + n (j + ]^,m\ . (5.47) 

Since Dunne and Halliday ^lO^ (working with coupling constant A = Sg) com- 
pute their characteristic polynomials Pm{E, A) in the subspace associated with 
decomposition (|5.45p , their polynomials Pm {E, A) factorize into our polynomials 

Rm{e,\) UM- 

J\//2 

[^Af(i?,A)]^^8^ = const W [R2,{E,g)r . (5.48) 

3 = {M/2} 
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For example, 

Pi = const -Ri , 
P2 = const R0R2 , 
P3 = const i?ii?3 , 
Pi = const RlRlRi . 

In this way the sl{2) algebra explains the "miracles" observed in [10]. 



6 Large orders of the perturbation theory 

6.1 Disappearance of the factorial divergence at negative 
even D 

In a general case the perturbative series for the energy of the ground state of 
the Z)-dimensional quartic anharmonic oscillator (jl.9p 



E{g,D)^J2^^'^(D)9'' 



(6.1) 



fc=0 



has factorially growing coefficients [7]: 



E^''\D) 



(-i)'^-+^r k 



D 



3^+^ 



2^/2 



6fc V3 2 4 



2)" tj:T{D/2) 

Oik-') 



(6.2) 



However, at negative even values of D function [F {D/2)] vanishes. Note that 
[F (D/2)]" appears on the RHS as a common factor so that at negative even 
values of D all terms of the systematic expansion in 1/fc become zero. This 
hints that at D = 0, —2, —4 the factorial growth of perturbative coefficients 
disappears. The same happens in for the potential 



Vir') = ir'+gr'^. 
In this case the factorial structure is different [1] 

1 (JV-1)^/^ 
n T{D/2) 



(6.3) 



^(fc)p)fc^oo [fc(7v_l)]! 



(fe) 
= (A) 



D/2 



X /fc(^/2)-l <^ 



1 J^ llv^ 



^r^(i^) 



Af-l- 



[l + 0(fc)] 



(6.4) 
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but we still have the same factor of [F {D/2)\ ~ . 

This factor of [r(D/2)]~ also appears in the asymptotic expression for 
E^'^^D) in the 0{D) symmetric model studied in Ref. [28]. This model de- 
scribes a particle on a D-dimensional sphere with a potential breaking the 
0{D + 1) symmetry of the sphere to 0{D). 

One can easily trace the origin of this universal factor [F (-D/2)] . In the 
semiclassical path integral derivation of large-order asymptotic formulas, the 
common factor of [F (£)/2)] comes from the integration over rotational col- 
lective coordinates of the classical solution breaking the 0{D) symmetry of the 
problem to 0(-D — 1). As a result of this integration, the final result is propor- 
tional to the surface area of {D — l)-dimensional sphere in the Z?-dimensional 
space: 

2 7r^/2 

^^ = T{Dj2) ■ (^-'^ 

Let us study the asymptotic behavior oi E^''\D) at negative even D = —2M. 
We will consider the case when the potential has the form 



V{r) 



—r 
2 



gU{r^) 



(6.6) 



where t/(r^) is some polynomial. We know that at I? = — 2M the energy of 
the ground state (understood as usual in the sense of analytical continuation) 
is given (at least for small g) by one of the roots of the algebraic equation 
(|1.3p . Since RM{E,g) is a polynomial in both g and E, the resulting depen- 
dence of E{g, —2M) is regular for most values of g. The radius of convergence 
of the perturbative expansion of E{g, —2M) is controlled by the singularity of 
E{g, —2M) closest to g = 0. These singularities appear at values of g corre- 
sponding to degenerate roots of polynomials RM{E,g) which are described by 
equations 

:i?A/(i?, go) 



-Rm(£^o,5o) 



dE 



==0. 



(6.7) 



E—Eq 



Note that a,t g — polynomial Rm{E, g) is given by expression (|2.42p which has 
no degenerate roots so that the point g ~ Q is regular and wc have a nonzero 
convergence radius. Equations ()6.7p may have several solutions. One has to 
choose the solution corresponding to the singularity closest to the point 5 = 
(on the main Riemann sheet). In the vicinity of the degenerate point goj Eq we 
can use the Taylor expansion 



RM{E,g)^UE-E, 



dE^ 



RMiE,go) 



(g-go) 



E—Eq 



d_ 
d'g 



Rm {Eq , g) 



J g=aa 
(6.8) 
Now equation (|1.3p leads to the following expression valid in the vicinity of the 
singularity of E{g) aX g — g^ is described by equation 



E{9) 



Eo- 




(6.9) 
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where 



c = -2.90 
Expanding 



—RM{Ea,g) 



g=ao 



dE^ 



RMiE,go) 



E — En 
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1 vr(fc-i)^g 



2^^T(k + l) V.go 



k=0 



and using the limit k ^ oo 



r(fc + i) 



fc-3/2 , 



we find from eq. (|6.9I 



EWiD)'=^l^l?^k-y'[g,m 



-fe 



(D = -2M). 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



6.2 Cases D = 0, -2, -4, -6 



Now we want to concentrate on the case of the quartic anharmonic oscillator 
p.9p and to consider cases D = 0, —2, —4, —6. 

The values D = and D = —2 are trivial because the corresponding poly- 
nomials i?o (|3.13p and i?i (|3.14p do not depend on g so that the energy is given 

by 



E{Q)^0, 

E{-2)^-l 

and almost all perturbative coefRcients E^'''> vanish: 

^('=)(0)=0 

£;W(-2) = -4o 



(6.14) 
(6.15) 



(6.16) 
(6.17) 



In the case D = —4 we insert polynomial R2{E,g) (|3.15p into eq. (|6.7p and 
find two solutions 



Eoi~4)^T^, .go(-4) = ±3-3/2 . 



(6.18) 



Taking into account that we are interested in the energy of the ground state 
that corresponds to 



D 

"2" 
i?(0,-4) = -2, 



[Eig,D)]^^, 



(6.19) 
(6.20) 



29 



one can check that the relevant singularity corresponds to the upper signs in eq. 
((6A81) 

i?o(-4) = --|, 5o(-4) =. 3-3/2 . (g21) 



Now we find from eq. (|6.10p 
and insert this into eq. (|6.13p 



c(-4) = - (6.22) 



E(k) (-4) ''^^lJlj2 k-^'^i^'^'^g'' (6.23) 

'^ fc>i 

for the coefficients E^^'' (D) of expansion (jl.Sp at Z? = —4. 

Other cases of negative even D can be analyzed in the same way. For exam- 
ple, solving equation (|6.7p for the polynomial R3{E,g) given by eq. (|3.16p we 
find for D = —6 

go(-6) = ^^=^, (6.24) 

36./3(5 + 2VT3) 



Eo{-6) = --^i (5 + 2V13) , (6.25) 



which results in the asymptotic behavior 



i,W(_6) '^^ 1^1 (^5 - :^)k-'/' iaoi-G)]-' ■ (6.27) 

6.3 Roots of polynomials E^^\D) 

Coefficients E'^'^-' {D) of the perturbative expansion (|6.fp are polynomials of de- 
gree k ^ 1. Let us denote the fc -f f roots of these polynomials v^^r'- 

e'^^'^ {vk,r) = (1 < r < fc + f ) (6.28) 

According to eqs. (|6.16p . (|6.f7p . at A: > 1 polynomials E^^^{D) have a common 
factor of £)(D-t-2), i.e. 

E^''\D) ^ D{D + 2)Pk{D) (6.29) 

where Pk{D) is a polynomial of degree fc — 1 with fc — 1 roots Vk^r- This means 
that the roots Vk^r include values and —2: 

{j^fc,r} = {0,-2,...} (6.30) 
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Note that the function [T (D/2)]"^ on the RHS of eq. ([^ has zeros at 

£> = 0, -2, -4, -6, . . . (6.31) 

The roots of the asymptotic expression (|6.2p at I? = 0, —2 appear exphcitly 
in the factorized expression on the RHS of eq. (|6.29p . The other roots should 
appear asyniptoticaUy at large k. In other words, in the full set of roots i'k,r 
there must be subsets converging to values D — —4, —6, . . . 
Let us study the first subset converging to 

D = -2M (M = -2, -3, . . .) 

We label this subset of roots with the index r — r{k, —2M): 

E^ {Vk,r(k-2M)) — 0, 

lim Vk,r{k.-~2M) = -2M . (6.32) 

k — ^oo 

We want to derive a large-fc asymptotic formula for 2M +Vkr(k,~2M) ■ To this 
aim we need an asymptotic formula for E'^^\D) in the double limit D -^ — 2M 
and fc — > oo. This asymptotic formula has the form 



ijWp) '=— ^-2Af (_i)/c+ir ^ ^^ 3fc+^^^^^(fe)(_2M). (6.33) 



2/ '-r(f) 

Note that we dropped the l/fc correction present in eq. (|6.2p but added the 
extra term E'^^\~2M) which becomes essential at _D ^ —2M because of the 
suppressing factor [F (-j)] i^ the leading term. At Z? ^ — 2Af we can simplify 
eq. Km 

E('^){D) ''^--■:^-2*^ (_i)fe+iA:!fc-M-i3^-*^^-y^ + £;(^-)(-2M) . (6.34) 

Here 

D\ D^-2M (-1)*' 1 



M 



r ^ ^^ yr^^TT, ■ (6-35) 

Inserting this and eq. (|6.13p into eq. (|6.34p . we obtain 

TT \ 2 

+ y'-^^k-^^'[goi~2M)]-\ (6.36) 

Using this asymptotic formula, we solve equation 

E'- ' {Vk,r{k-2M)) = (6.37) 
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and find 



■,M 



fe— (-6) / , o.n (-l) 



^.,.(fc.-2M) + 2M '^ i^v/'rc(-2M)^fc^'^-(V2) [35o(-2Af)] 



-fc 



(6.38) 

Using above expressions (O^ . ((O^ . (I^?^ . (p?^ . for c(-2A/) and go{--2M) 
at M = 2, 3, we obtain 

^fe,.(fe,-4) + 4 '^ 12v^H)!a:3/23'=/2 , (6.39) 

5\/l3- 1 



i^k,rik,-6) + 6 "^^ 12 J27r ( 5 - -^) ^^fc5/2 



12J3 (5 + 2^/T3) 



\/T3; fc! 

(6.40) 
Asymptotic formulas (|6.39p . (|6.40p agree with results of the direct numerical 
calculation using methods of Refs. [H |2l [5l [7] . The roots approaching —4 appear 
rather early. They are listed in Table [1] For small odd values k — 5,7,9 these 
roots have small imaginary parts. Starting from k = 10, the roots Vk,r(k,i) close 
to —4 become real and exhibit a very fast factorial convergence to —4 which 
is described by the asymptotic formula (|6.39p . Already at fc = 11 this formula 
works with accuracy better than 7%. 

At large values of k there appear roots converging to other integer values 
D — —6, —8, . . . The roots approaching the value —6 are listed in Table [2l At 
large k they are well described by asymptotic formula ()6.40p . 

6.4 Distribution of roots in the complex plane 

One should keep in mind that apart from the stable roots at Z? = 0, —2 and the 
roots asymptotically approaching points D = —4, —6, . . ., there are many other 
roots i'k,r- Most of them are complex. In order to get an impression about 
the general distribution of roots in the complex plane, we show them for the 
cases fc = 7 (Fig. [3l), fc = 8 (Fig. [l, and fc = 30 (Fig. E]) and fc = 60 (Fig. 
[S]). From the plots for fc = 30 and fc = 60, it is clearly seen that at large fc 
there appears a certain regular structure in the distribution of complex roots. 
It is an interesting problem to give a complete asymptotic description of roots 
in the complex plane. Here we make only preliminary comments about the 
distribution of roots. 

It is convenient to work with polynomials Pk{D) (|6.29p . Using eq. (jl.lOp . 
we find: 

Pk{D) "^^ (-l)'=+ir ( fc + - I 3"+^^-- . (6.41) 

^ ' ^ ^ V 2/ 7rr(i:>/2 + 2) ^ ' 

The polynomial Pk{D) has roots iyk,r which coincide with the roots of E'^'^\D) 
up to the two stable roots 0, —2: 

{vks,Vk,2^- --^k^k+i} = {vk,i,i'kai- ..I'fe.fe-ijO, -2}. 
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k 


Vk,r{k.-i) + 4 


k 


Vk.r{k-i) +4 


5 


-3.22834 ±i0.426293 


33 


-4.79523448 x 10"^'' 


6 


-3.44545 


34 


+2.55611474 X 10"^^ 


7 


-3.63083 ±i0.34226 


35 


-1.32186040 X 10-^» 


8 


-3.76443 


36 


+6.63762841 x IQ-^'J 


9 


-3.9583 ±i0.226557 


37 


-3.23912410 X lO^^i 


10 


+0.04231592827 


38 


+ 1.53735607 X 10"^^ 


11 


-0.01231265412 


39 


-7.10198670 X 10"^'' 


12 


+0.00178433080 


40 


+3.19560318 X 10"^^ 


13 


-0.00027422590 


41 


-1.40147930 X 10-^'' 


14 


+0.00003787462 


42 


+5.99459657 x lO"^" 


15 


-4.86252995 x lO"" 


43 


-2.50229008 x 10"^'^ 


16 


+5.80950053 X 10" "" 


44 


+ 1.01993415 X 10"""' 


17 


-6.49387664 x 10"** 


45 


-4.06166432 x 10"^^ 


18 


+6.81906230 x 10"^ 


46 


+ 1.58111340 X 10-''^ 


19 


-6.75145346 X lO^^" 


47 


-6.01961315 X 10-4^ 


20 


+6.32321355 X 10"^^ 


48 


+2.24249068 x lO"*" 


21 


-5.61842521 x lO"^"" 


49 


-8.17805838 x IQ-** 


22 


+4.74864227 X 10"^^ 


50 


+2.92092144 x IQ-^^ 


23 


-3.82680164 X lO"^"* 


51 


-1.02217301 X 10-'^" 


24 


+2.94682238 x lO^^'^ 


52 


+3.50623616 x lO"'^^ 


25 


-2.17261017 X 10-1^ 


53 


-1.17934411 X 10-'^-^ 


26 


+1.53640755 x 10"" 


54 


+3.89122545 x lO"'^'^ 


27 


-1.04388321 X lO^^** 


55 


-1.25990087 X 10"^" 


28 


+6.82474923 x 10"^" 


56 


+4.00444478 x 10"'^** 


29 


-4.29961693 x lO^^^ 


57 


-1.24982870 X lO"'^*^ 


30 


+2.61370107 X 10-^^ 


58 


+3.83179265 x lO^^i 


31 


-1.53497157 X 10"^^ 


59 


-1.15433788 X 10"^^ 


32 


+8.71891707 X 10"^'^ 


60 


+3.41802128 X lO-*^" 



Table 1: Roots I'k.rikj-A) of polynomials Pd{N) approaching the value —4. 
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k 


'^k,r(k.-6) + 6 


k 


'^k,r{k.-6) +6 


21 


+0.03432898 


41 


+1.44684 X 10-i« 


22 


-0.011016 


42 


-1.55969 X lO-i'^ 


23 


+0.0020276 


43 


+1.63996 X 10-^" 


24 


-0.00040899 


44 


-1.68294 X 10-"i 


25 


+0.000076759 


45 


+1.68654 X 10-^^ 


26 


-0.000013876 


46 


-1.65141 X lO-""-^ 


27 


+2.40442 X 10-^ 


47 


+1.58077 X lO-^'^ 


28 


-4.00519 X 10"' 


48 


-1.47999 X lO-""^ 


29 


+6.42173 X 10-** 


49 


+1.35591 X 10-"" 


30 


-9.92472 X 10--^ 


50 


-1.21614 X 10-"' 


31 


+1.48040 X 10--^ 


51 


+1.06834 X 10-"** 


32 


-2.13385 X 10-^" 


52 


-9.19592 X 10-^" 


33 


+2.97550 X 10-^^ 


53 


+7.75910 X 10-^1 


34 


-4.01814 X 10-^^ 


54 


-6.41992 X 10-^" 


35 


+5.26006 X 10-^^ 


55 


+5.21090 X 10-^^ 


36 


-6.68131 X lO-i'^ 


56 


-4.15065 X lO-^'^ 


37 


+8.24173 X 10-1^ 


57 


+3.24558 X 10-^^ 


38 


-9.88147 X 10-1'' 


58 


-2.49222 X 10-^" 


39 


+1.15242 X 10-1'' 


59 


+1.87991 X 10--^^ 


40 


-1.30830 X 10-1' 


60 


-1.39342 X 10-^« 



Table 2: Roots Vk.r(k-&) of polynomials Pk{D) approaching the value —6. 
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We can represent Pk{D) in the form 



fc-i 



Pk{D) = hX{{D - Vk.r 



r=l 



The coefficient /3k was computed in Ref. [6]: 

p/ 3fc-l -> 

f3k = (-i)'=+i2'=-2 — ^y 2 ) 



(fc + i)!r(^) 



We find from eqs. ^Z3l, (p^ 

fc-i 



PtID) 



r— 1 ^ ' 



m 



D/2 



(6.42) 



(6.43) 



(6.44) 



PkW ^J^\ h,rj TiD/2 + 2) 

Differentiating the logarithm of this expression with respect to D, we obtain 
fc-i 



E 



1 fc^oo 1 

D - Vk,r 2 



In {Qk) - V ( — + 2 



where V'(-2^) is the logarithmic derivative of the F function 

r'(z) _ '^ / 1 1 



r(2;) 
and 7 is Euler's constant. Thus 



fc-i 

n 



1 



1 



-D - i>fc, 



In (6/s) 



E 

fc=0 



7 + E 

fe=0 



z + k 



D + A + 2k 



1 



(6.45) 



(6.46) 



(6.47) 



This relation requires the appearance of roots Vk.r converging to —4, —6,. . . but 
does not forbid the existence of additional complex roots going to infinity or 
making a quasicontinuous distribution in the limit k —^ oo. Numerical calcula- 
tions show a large amount of complex roots at large k (see Figs. [5]and[n|). 

In the limit of large k, one can compute the product of all roots i'k,r of the 
polynomials Pk{D) setting £> = in cqs. (|6.4ip . (|6.42p and using expression 
((05)) for f3k: 



k-l 



/3fe4^' 



Jl j.,^^ ^=^ L^T (k) "^^ i-l)''+'3k' 



\k+lou2 



;V3 



(6.48) 



This gives us an asymptotic estimate for the geometric average of all nonzero k 
roots i'k,r of Ek{D) (including the root Vk^r = —2) 



iWA) 



Vk.^T^O 



■ /fc-1 

2 n i^fc.> 



\r=l 



l/fc 



'--oo k ,^.1/k 



3n/3 



(6.49) 



The growth of this quantity with k shows that most of the roots go to infinity. 
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Figure 3: 8 roots v^^r of the polynomial E'^'^\D). Two complex conjugate roots 
—3.63 ± iO.34 are in the vicinity of the negative even value D = —4. 
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Figure 4: 9 roots i^g^r of the polynomial E'-^^{D). One root at —3.76443 is close 
to the negative even value —4. 
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Figure 5: 31 roots t^so.r of the polynomial E^^^^D). One can see a discrete set 
of roots close to even negative values up to D = —8. Most of the roots belong 
to the quasicontinuous set formed in the complex plane. 
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Figure 6: 61 roots feo.r of the polynomial E^^^'>{D). One can see a discrete 
set of roots close to even negative values up to Z? = —12. There is also a 
quasicontinuous set formed in the complex plane. 
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7 Conclusions 

Thus we have estabhshed a close connection between 

1) exactly solvable features of the anharmonic oscillator in negative even 
dimensions, 

2) disappearance of the factorial growth of perturbative coefficients in neg- 
ative even dimensions, 

3) fast inverse factorial convergence of roots of polynomials E^'''> (D) to the 
negative even points. 

Although the content of this paper was restricted to quantum mechanics, the 
main motivation comes from quantum field theory. In contrast to quantum me- 
chanics where the large-order behavior of the perturbation theory can be easily 
tested using direct numerical calculations in cases when the analytical methods 
fail or raise doubts, in quantum field theory the power of both analytical and 
numerical tools is rather limited. The slow 0{k~^) convergence of perturbative 
coefficients to the asymptotic form [see e.g. eq. (|6.2p ] is rather disturbing in 
quantum mechanics but one still can reach the asymptotic regime in large or- 
ders. In quantum field theory calculations rarely go beyond four or five loops. 
In this situation the construction of quantities whose perturbative expansion 
has a fast factorial convergence to the asymptotic form is very important. 

Polynomials i?^*^^ (D) and their roots have analogs in quantum field theory. 
In many field theoretical models we have dependence on various parameters 
(e.g. number of colors, flavors etc.) and Feynman diagrams in any order lead 
to a polynomial (or fractional polynonnal) dependence on these parameters. 
The control of this polynomial dependence is usually trivial compared to the 
hard work needed for the calculation of loop integrals. Therefore available 
multiloop results in quantum field theory provide many opportunities for testing 
the roots of these polynomials. There is some evidence that the inverse factorial 
convergence of roots found in the case of the anharmonic oscillator has analogs 
in some field theoretical models, e.g. in the A^-component 0{N) symmetric 
model [27] . An indirect theoretical argument comes from the presence of an 
inverse F function in field theoretical asymptotic expressions. This inverse F 
function may generate zeros similar to the zeros at D = 0, —2, —4 . . . coming 
from the factor [F(Z?)]~ in eq. (|1.10p in the case of the anharmonic oscillator. 
Since the convergence is very fast, it may happen that one can detect it (or its 
signals) in available multi-loop perturbative expressions. 

Acknowledgments. I appreciate discussions with G.V. Dunne, L.N. Lipa- 
tov, V.Yu. Petrov and A. Turbiner. 
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